Quenching along a gapless line: A different exponent for defect density 
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We use a new quenching scheme to study the dynamics of a one-dimensional anisotropic XY spin- 
1/2 chain in the presence of a transverse field which alternates between the values h -\- 5 and h — 5 
from site to site. In this quenching scheme, the parameter denoting the anisotropy of interaction 
(7) is linearly quenched from —00 to +00 as 7 = t/r, keeping the total strength of interaction J 
fixed. The system traverses through a gapless phase when 7 is quenched along the critical surface 
h? = 5^ + in the parameter space spanned by h, 5 and 7. By mapping to an equivalent two- 
level Landau-Zener problem, we show that the defect density in the final state scales as a 
behavior that has not been observed in previous studies of quenching through a gapless phase. We 
also generalize the model incorporating additional alternations in the anisotropy or in the strength 
of the interaction, and derive an identical result under a similar quenching. Based on the above 
results, we propose a general scaling of the defect density with the quenching rate r for quenching 
along a gapless critical line. 
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I. INTRODUCTION 

The dynamics of a quantum system swept across a 
quantum critical point at a uniform rate has been stud- 
ied extensively in recent years. Since a quantum phase 
transition^'^ is necessarily accompanied by a diverging 
correlation length as well as a diverging relaxation time, 
the dynamics of the system cannot be adiabatic for the 
entire period of the evolution however slow the varia- 
tion in the parameter may be. (The relaxation time of 
a quantum system is given by the inverse of the energy 

\ gap which goes to zero at the quantum critical point). 
Assuming that the system was initially prepared in its 

' ground state, the non-adiabaticity near a quantum criti- 
cal point prevents the system from following its instanta- 
neous ground state resulting in the production of defects 
in the final state. 

The Kibble-Zurek (KZ) argument^ asserts that the 
non-adiabatic effect becomes prominent only close to the 
critical point when the rate of change of the Hamiltonian 
is of the order of the relaxation time of the underlying 
quantum system. When a parameter of the quantum 
Hamiltonian is varied as t/r, where r is the character- 
istic time scale of the quenching, the above argument 
predicts a density of of defects in the final state that 

1 scales as \ J t'^^ / ^'^'^+'^) in the limit of r ^ 00. Here v 
and z denote the correlation length and dynamical expo- 
nents, respectively, characterizing the associated quan- 
tum phase transition of the d-dimensional quantum sys- 
tem. The KZ prediction has been verified for various 
exactly solvable spin models when quenched across a 
critical^ or a multicritical® point at a uniform linear 
rate. The above studies have been generalized to explore 
the defect production in a non-linear quench across a 
quantum critical point where a parameter in the Hamil- 
tonian is varied as It/rj" with a > 0^. Recent experimen- 
tal studies on the dynamics of quantum systems^, espe- 



cially quantum magnets^, ultracold atoms trapped in op- 
tical latticesi^ and spin-one Bose-Einstein condensates^, 
have paved the way for a plethora of related theoretical 
studiesi^i^i^i^^i^^. 



Another interesting scenario emerges when a low- 
dimensional quantum system is quenched through a gap- 
less phase or an extended quantum critical regioni^ii^i^S. 
It has been established that when a d-dimensional system 
is quenched along a (d — m)-dimensional critical surface, 
the scaling of the defect density with r is modified to a 
generalized KZ form given by i/-7-'"''/(^f+i)i9. 



In the present work, we explore the dynamics of a 
one-dimensional anisotropic XY spin-1/2 chain in the 
presence of a transverse field which alternates between 
h -\- 5 and h — 5 from site to site. We employ a new 
quenching scheme in which the parameter determining 
the anisotropy of interaction is quenched as t/r, keeping 
the strength of interaction fixed, in such a way that the 
system is driven along a gapless line on a critical surface 
in the parameter space. We show that the density of de- 
fect scales as a result that cannot be explained 
by the previous studies on the quenching through a gap- 
less phase. We also propose a general scaling relation for 
such a quenching dynamics along a gapless line. 



The paper is organized as follows. Our model, the 
quenching scheme and the results obtained for the gen- 
eration of defects are presented in Sec. II. At the end of 
that section, we propose a general scaling relation for the 
defect density when a system is quenched along a gapless 
line. We end with some concluding remarks in Sec. III. 
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II. THE QUENCHING DYNAMICS AND THE 
RESULT 



The Hamiltonian of the spin- 1/2 anisotropic XY 
model with an alternating transverse field is given 



1=0.5 J=l 



H 



+ (J, - Jy){a^a^+, - a]a)^^) + {h - (-lF<5)a|)], (1) 

where cr's denote the Pauli spin operators satisfying the 
standard commutation relations. The strength of the 
transverse field coupled to the operator cr^ alternates be- 
tween h + S and h — S on the odd and even sites re- 
spectively. We have chosen all the interactions and the 
fields to be non-random. Henceforth, we shall refer to 



Jx + Jy — J SiS the strength and j — Jx 



Jy as the 



anisotropy of the nearest-neighbor interaction. 

To map the spin operators to spinless fermion oper- 
ators using the Jordan- Wigner transformatio] 



.23.24 



we 



note that the presence of two underlying sub-lattices ne- 
cessitates the introduction of a pair of fermion operators 
a and for even and odd sites as defined below 



'2j 



'2j+l 



J-1 

62 exp[z7r 



1=1 
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=1 



a;, Li exp[z7rE4/^2/ 



J-1 
1=0 



ia2i+i]- (2) 



Using a restricted zone scheme (where the wave vector 
k ranges from — 7r/2 to 7r/2) in the Fourier space, the 
Hamiltonian can be written as 



H = j2h, = 



k 



where Ak is the column {a\,a-k^b\,b^k)- The 4x4 
Hermitian matrix is given by 

h + Jcosk i"fsmk —6 

—27 sin /c —h — Jcosk S 
6 J cos k — h 17 sin k 

—5 — i7sinfc —Jcosk + h 



(3) 



The excitation spectrum of the Hamiltonian H is now ob- 
tained by diagonalizing the reduced Hamiltonian matrix 
Hk and is given by 



±2 



cos^ k + j'^ sin^ k 



s^k + d^j'^ sin^ /c]i/2. 



(4) 



Denoting the four eigenvalues by ±A^, we can write the 
spectrum of the Hamiltonian in the form 



-7r/2<fc<7r/2 



(5) 



'^\\S"= h 
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FIG. 1: Phase diagram of the XY chain in an alternating 
transverse field. We have chosen J = 1 and the critical lines 
are drawn in the 7 = 0.5 plane. Two special points (h = ±J, 
(5 = 0) and {h = 0,5 — ±7) are shown on the phase bound- 
aries. The spin chain undergoes a quantum phase transition 
with v = 2 and z = 1 when these points are approached along 
the dashed line. On the other hand, the dotted line shows the 
direction of the quenching of the transverse field. We quench 
the system along a gapless line parallel to the 7-axis (perpen- 
dicular to the plane of the paper) passing through the phase 
boundary = + J^. 



where ^ is the quasiparticle creation operator corre- 
sponding to the mode (fc, j/). In the ground state the 
levels — and — are filled. At the quantum critical 
point, A^ vanishes at an ordering wave vector, and the 
critical exponents are obtained by studying the behavior 
of AjT in the vicinity of the critical point. 

The minimum energy gap in the excitation spectrum 
occurs at fc = and k = 7r/2. The corresponding phase 
boundaries given by h"^ — 5^ + ,P and (5^ = ft-^ -f 7^ 
signal quantum phase transitions from a paramagnetic 
to a ferromagnetic phase and a dimer to ferromagnetic 
phase respectively (see Fig. 1). Let us define a new set 
of Pauli matrices r as 



and 



(-l)Vf 



so that the commutation relations of the Pauli matrices 
are preserved. It is interesting to note that under this 
unitary transformation, we arrive at a set of duality re- 
lations given by /i — !■ —(5, 5 —h, J —7 and 7 — > — J; 
this signifies that the ferro-para transition and the ferro- 
dimer transition at ^ + and 5^ — h'^ + "/^ , 
respectively, are essentially identical, both belonging to 
the quantum Ising universality class^''' with ly = z = 1. 
The phase boundary given by h'^ = + with 7 ar- 
bitrary and J held fixed, defines a critical surface in the 
parameter space spanned by /i, <5 and 7. Similarly, the 
phase boundary — /i^ +7^ with arbitrary J once again 
defines another critical surface when 7 is held fixed. For 
7 = 0, a gapless phase exists with an ordering wave vec- 
tor cos A: = VW^~P/J for < h"^ < 5"^ + J"^ . The sys- 
tem undergoes a quantum phase transition from a gapless 
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phase to a gapped phase when \h\is increased beyond the 
critical value given by he = \/S^+J^. 

The special case with S — refers to the well studied 
anisotropic XY spin- 1/2 chain in a transverse field^. In 
this model, there exists an Ising transition line at h = ± J 
from the ferromagnetically ordered phase to a quantum 
paramagnetic phase. On the other hand, the vanishing of 
the gap at 7 = signifies another quantum phase tran- 
sition belonging to a different universality class between 
two ferromagnetically ordered phases. 

It can be established using a numerical diagonalization 
of the time-dependent Schrodinger equation involving the 
reduced Hamiltonian matrix that when the transverse 
field h or the alternating term S is quenched as t/r from 
— oo to oo, so that the system crosses the quantum critical 
lines as shown in Fig. 1, the density of defects n in the 
final state satisfies the Kibble-Zurek predictior^^^ . Our 
interest however lies in the generation of defects when 
the system is quenched along a gapless line. To achieve 
such a quenching, we vary the anisotropy parameter 7 
linearly as 7 = t/r from —00 to 00, keeping h, S and J 
fixed in such a way that the system always lies on the 
phase boundary — + J^. In the limit t — s- —00, 
7 is large and negative and hence in the ground state. 



the expectation value 



-1. On the 



other hand, for an adiabatic evolution during the entire 
period of dynamics, this expectation value should be -1-1 
in the final state. One can choose different critical lines 
for 7-quenching by choosing different values of h and 5 
on the critical surface. 

The eigenvalue given by Eq. ^ can be written as 



A- = [{h - yJS^ + ,P cos2 kf + 2hyJP + ,P cos2 k 



+ 7^ sin^ k - 2^h^S^ h^,P cos2 k + S^i^ sin^ fej^/^. 

(6) 

On the gapless line /i^ = -|- J^, the dispersion of the 
low-energy excitations at fc — s- can be approximated as 



J4fc4 



7^ J2fc2 



V 4(52 + J2) ' 52 + J2 



(7) 



A close inspection of the above excitation spectrum 
suggests that when 7 is quenched along the gapless line, 
Eq. d?]) can be mapped to the spectrum of a 2 x 2 Landau- 
Zener (LZ) Hamiltonian^^i^l with two linearly approach- 
ing time-dependent levels. To show this explicitly, we 
note that in the limit of very slow quenching, r — > 00, 
defects are produced by sets of modes between which the 
energy gap is very small. For the Hamiltonian in Eq. l[3|. 
this occurs in the region of fc = if we take \\} = (5^ _)_ 
Let us first set 7 = and fc = 0. We then see that there 
are two modes, called |/) and |//), whose energies are 
zero; the other two modes have energies which are both 
far from zero and far from each other, and can therefore 
be ignored in a slow quenching calculation. The zero 



energy modes are given by 

I/) 



+ + J)2 



\n) = 



5 



h + J, 



5 

JS^ + ih + J)2 I h + J 
* 



1 



(8) 



We now deviate slightly from fc = 0, still keeping 7 = 
in Eq. ^ . Doing degenerate perturbation theory to first 
order in Jfc^, we find that the modes |/) and \II) remain 
eigenstates of the Hamiltonian, but their energies are now 
given by 



21.2 



E+ ^ ± 



J'k 



(9) 



where we have used the relation h? = 5"^ + ,P . Finally, 
we introduce the terms involving 7 in Eq. ([3]). To first 
order in 7, we find that the Hamiltonian in the basis of 
|/) and \II) is given by 



1 



ijJk 



—i"fjk 



(10) 



If we now perform a unitary transformation and vary 7 
in time, we see that the Hamiltonian is of the LZ form. 



hk = 



^{t)k j2fc2/2 
j2fc2/2 -7(t)fc 



(11) 



where 7 and J are renormalized parameters given by 
7 = -iJ/^fPTl^ and J2 = J2/V(52 + J2. The diago- 
nal terms in Eq. pT|) describe two time-dependent levels 
approaching each other linearly in time (since 7 = ^/r), 
while the minimum gap is given by the off-diagonal term 
j2fc2/2. The probability of excitations pk from the 
ground state to the excited state for the fc-th mode is 
given by the Landau-Zener transition formul a^^i^'^ 



Pk 



exp 



27rJ'*fc 



8fc 



dt 



exp 



nJ^k^T 



(12) 



Note that for large t, pk is dominated by values of fc close 
to 0. The density of excitations n in the final state and 
in the large r limit is obtained by integrating over all 
modes in Eq. P^. 



2 

TT 



n/2 



dk Pk 



1 



-1/3 ■ 



(13) 



The numerical integration of Eq. for 6 = J = 1 is 
shown in Fig. 2 (a) ; although we have used the expression 
given in Eq. for all values of fc, the dominant contri- 
bution to n comes from the region near fc = where (jl2[) 
can be trusted. This shows that when quenched along a 
gapless line, the density of defects in the final state ex- 
hibits a slower decay with t as compared to the 1/\/t 
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followed by a Fourier transformation, we find an excita- 
tion spectrum of the form 

A± = [h^ + S^ + (J2 + 7^) cos^ fc + 72 sin^ k 

- 2^{h^ + 72 cos2 k){6^ + J2 cos2 k) + (5272 sin^ k]^/^, 

(14) 

where the eigenvalue has to be analyzed to explore 
the quenching dynamics. Eq. (fT4|) shows that the role 
of the alternation 7^ is to renormalize the strength J; 
consequently, the phase boundary separating the para- 
magnetic and the ferromagnetic phase gets shifted to 
/i2 = 52 -I- j2 — 72 ^ with arbitrary 7, at which the gapless 
excitations occur ai k — 0. Quenching 7 as t/r along 



the new phase boundary h 



7, 



2 with fixed 



values of J, h and 7^, once again takes the system along 
a gapless line on a critical surface. 

The dynamics can be reduced to a two-level problem as 
before, and the low-energy excitations above the gapless 
line are given by 



FIG. 2: Fig. 2 (a) shows the variation of defect density n 
with T for 5 — J — 1 obtained by numerically integrating Eq. 
(|13[) . using pk given in Eq. (12). Fig. 2 (b) shows n obtained 
by direct numerical integration of the Schrodinger Equation 
with 5 — 0. In both figures, the fitted data shows a slope of 
— 1/3 confirming the n ~ behavior. 



which is observed in the case when the gapless line is 
crossed by varying hS^. Our study also reveals that al- 
though the spectrum of this system involves four levels, 
the quenching dynamics along the gapless line is essen- 
tially a two-level problem which can be studied using an 
effective Landau-Zener theory with parameters renormal- 
ized by S. 

The case 6 = corresponds to quenching the system 
along the Ising critical line of the transverse XY model, 
and one obtains an identical scaling of the defect density. 
Fig. 2 (b) shows the 1 /3 power-law obtained by numeri- 
cally solving the Schrodinger equation for the anisotropic 
transverse XY model when the anisotropy parameter 7 
is quenched along the gapless line h = —J = 1. (The 
Hamiltonian in Eq. ^ has a 2 x 2 block diagonal form 
if (5 = 0). One can also propose an alternative quenching 
scheme where the strength J is quenched from — cxd to 
00, keeping h, S and 7 constant with S"^ = + j^. The 
duality relation discussed above leads to the conclusion 
that this quenching scheme is equivalent to the previous 
one with 7 = t/r, and it yields a similar I/t^^^ behavior. 

The XY chain with an alternating transverse field can 
be further generalized by incorporating additional alter- 
nations in the strength or in the anisotropy of the interac- 
tion with the period of alternation being two. We denote 
the alternation in the strength and the anisotropy by Js 
and 7s, respectively, and for simplicity choose = 0. 
Using a similar Jordan- Wigner transformation (Eq. Q) 



'(J2-72)2/c4 72J2fc2 



4(^2 + J2) 



52 + J2 



(15) 



The defect density n in the final state decreases with r as 
as can be derived from the Landau-Zener formula 
in an identical fashion. Similarly, one may set 7^ = 
and Js ^ and consider an equivalent quenching scheme 
resulting in an identical scaling of the defect density. 

The behavior of the defect density when quenched 
along a gapless line suggests the following general scaling 
relation of the defect density for a d-dimensional quan- 
tum system. Let the excitations on the gapless quantum 
critical line be of the form ~ a|fc|^, where z is the dy- 
namical exponent and the parameter a = t/r is quenched 
from —00 to CXD. Using a perturbative method involving 
the Fermi Golden rule along with the fact that the sys- 
tem is initially prepared in the ground state, the defect 
density can be approximated as^^ 



d'^k 
(2^ 



da {k\-^\0) e^^^°'^"£("')'*"'|2.(16) 
da 



Assuming a general scaling form of the instantaneous ex- 
citation 6Ljj:{a') = k°-f{^), where k = and k"- de- 
notes the higher order term in the excitation spectrum 
on the gapless line. Defining a new variable ^ = ak^~°', 
we obtain the scaling behavior of the defect density as 



1/r 



d/{2a-z) 



(17) 



The case d = 1, a = 2 and z = 1 has been discussed 
in the present work. Note that the correlation length 
exponent v does not appear in the expression in Eq. (jl7p 
because our quench dynamics always keeps the system on 
a critical line. 
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III. CONCLUSION 

In conclusion, we have studied the defect density pro- 
duced in the final state when a generalized spin-1/2 XY 
chain with an alternating transverse field as well as an al- 
ternating nearest-neighbor interaction is quenched along 
the Ising critical line by varying the anisotropy parameter 
7. We show that the non-adiabatic transition probabil- 
ity and hence the defect density can be estimated using 
an equivalent Landau-Zener problem in which the pa- 
rameters are renormalized by the alternating parameters 
S and 7s (or Jg). We find that the defect density de- 
cays with the characteristic time scale of quenching given 
by The defect scaling exponent obtained here 

does not fit the KZ scaling 1/t'^'^/'-^'^^^\ In the present 
quenching scheme, the system is always on a critical sur- 



face, and therefore the critical exponent v does not ap- 
pear in the scaling of the defect density. The quenching 
scheme used here is different from the other quenching 
schemes through a gapless phas o"'^^'^^i^° where the system 
starts from a non-critical (gapped) point, goes through a 
critical (gapless) point or critical surface, and eventually 
ends again at a non-critical point. 
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